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We calculate numerically the sizes S of jumps (avalanches) between successively pinned config- 
urations of an elastic line {d = 1) or interface (d = 2), pulled by a spring of (small) strength 
in a random-field landscape. We obtain strong evidence that the size distribution, away from the 

small-scale cutoff, takes the form P(5) = f p(S/^.) where -.^ ^ ^ m'-C the scale of 
avalanches, and the roughness exponent at the depinning transition. Measurement of the scal- 
ing function /(s) :— s^p(s) is compared with the predictions from a recent Functional RG (FRG) 
calculation, both at mean-field and one-loop level. The avalanche-size exponent r is found in good 
agreement with the conjecture r = 2 — 2/(d-|-^), recently confirmed to one loop via the FRG. The 
function f{s) exhibits a shoulder and a stretched exponential decay at large s, ln/(s) ~ — s*, with 
5 « 7/6 in d = 1. The function /(s), universal ratios of moments, and the generating function 
(e"^") are found in excellent agreement with the one-loop FRG predictions. The distribution of local 
avalanche sizes S^, i.e. of the jumps of a subspace of the manifold of dimension d^, is also computed 
and compared to our FRG predictions, and to the conjecture = 2 — 2/[d^ + (). 



I. INTRODUCTION 

Elastic objects pinned by a random substrate are 
ubiquitous in nature. The competition between elastic 
restoring forces and quenched disorder results in multi- 
ple metastable states. Upon applying an external force 
one observes collective jerky motion which proceeds by 
sudden jumps, called avalanches. Examples are the 
Barkhausen noise in magnets yj [H [31 SI [51 IHl |7] , jumps in 
the creep motion of magnetic domain walls [HI [HI UHl IH] , 
avalanches in the depinning of a contact-line of a fluid 
[121 ESI [HI [ll]i or in dislocation and crack propaga- 
tion [ini HZl [iSl [ini, and stlck-sHp motion of e.g. tec- 
tonic plates, responsible for earthquakes [^ \Tn I23j. 
Avalanches have also been studied in models without 
quenched substrate disorder, such as in sandpile mod- 
els and in granular matter [Ml US [551 [27] . An impor- 
tant characteristics of avalanche motion is its scale in- 
variance, self-organized criticality, and a broad distribu- 
tion P{S) ^ S^"^ of the sizes S of avalanches, for sizes S 
between a small- and large-scale cutoff S'min <C 5 <C Sm. 
Pinned elastic manifolds are an important prototype of 
a much wider class of phenomena, reaching far outside 
physics, e.g. into economy and finance, where extreme 
(and sometimes catastrophic) events are sufficiently fre- 
quent and large to dominate most observables. In this 
context, it is clearly of importance to understand how the 
avalanche-size probability is cut off at the large scales, for 

Although avalanche motion of pinned manifolds has 
been studied for a while in numerics [28^ .29^ 30}, most 
work focused on measuring the avalanche-size exponent 
T, with minimal guidance from the theory. This is mainly 
because no analytic approach was available besides mean- 
field and scaling arguments. The most notable one was 
proposed by Narayan and Fisher (NF) [3T] on the basis 



of the unproved assumption that the avalanche density 
remains finite at the depinning threshold, resulting into 



Here is the roughness exponent at the depinning tran- 
sition. Progress both in constructing the field theory of 
the depinning transition [32l [33l [SH [355 following the pi- 
oneering work on the Functional RG (FRG) [3T1 [351 [37] 
and in developing new powerful algorithms [THl [35} [321 
l40l l4n l42l 143] had focused mostly on structural prop- 
erties of the pinned manifold, such as the precise deter- 
mination of C. Even an appropriate definition of static 
and dynamic avalanches, allowing contact with the field 
theory, had remained unclear. It was given in the stat- 
ics [351 SS] and at depinning [JT] [35] using a confin- 
ing quadratic potential; it led to the measurement, with 
great accuracy, of the renormalized disorder correlator 
A(u), i.e. the fixed point of the FRG. Only very recently 
we succeeded in computing the distribution of avalanche 
sizes within the FRG [IHl [SHI [5T|. The calculation at 
tree level gave mean-field predictions (some of them new 
and non-trivial) , valid above the upper critical dimension 
due — 4, i.e. for d > 4. The one-loop calculation gave an 
expansion to order 0(e), with e = 4 — d. Remarkably, the 
conjecture ([T| was confirmed to 0{e) accuracy. It is thus 
of great interest to test these predictions in numerics. 

The aim of the present paper is to compute numerically 
the jumps (avalanches) between successively pinned con- 
figurations of an elastic line {d = 1) and interface {d = 2). 
The convenient setting to compare with the recent pre- 
dictions from the FRG is to submit the manifold to an ex- 
ternal quadratic well, i.e. a spring. We will study mostly 
random-field disorder, but we also check that the results 
are the same for random-bond disorder, as is predicted at 
depinning and was checked in our previous work [48] for 



2 



the renormalized disorder correlator A(u). Most of the 
numerical method is similar to our previous work [48j. 

The outline of this article is as follows: We define in 
section JITl the model and numerical procedure; and in 
section IIIII an avalanche, its size, the characteristic scales 



and the scaling functions. In section IV the reader will 
find our numerical results for the avalanche-size distri- 
bution, and their comparison to our analytical results in 
d = 1 and d = 2. In section |V] we compare our numeri- 
cal and analytical results for the universal ratios of alge- 
braic moments, the r„. In section [VT] we do the same for 
the generating function of exponential moments, i.e. the 
characteristic function of the size distribution, den oted 
Z{X), in = 1 and d = 2. Finally, in Section (VII I, wc 



compute the distribution of local avalanches in d 
compare with the predictions. 
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II. NUMERICAL PROCEDURE: PARABOLA 
AND METASTABLE STATES 

Let us now describe the model and algorithm, in d = 1 
for simplicity. The procedure is very similar to our pre- 
vious work [48J. The interface is discretized as u(x) = Ui, 
i = 1, . . . , L, and periodic boundary conditions are taken: 
ua — ul, — ui. We start from a flat interface 

(ui = 0) embedded in a parabolic potential. The equa- 
tion of motion is 

dtUi = m^{w - Ui) + Ui+i + Ui^i - 2ui + F{i, u,) (2) 

F(i, Ui) is the disorder force. We distinguish two different 
microscopic disorders: 

(i) random force (RF): for each integer value of Ui we 
take a random number extracted from a normal distri- 
bution. The value of the random force for non-integer 
values of Ui is given by the linear interpolation of the 
forces at the two closest integers Ui. Forces for different 
i are independent. 

(ii) random bond (RB): the random force is derived 
from a random potential: F{i,Ui) = —duiV{i,Ui). For 
each integer value of Ui, the potential is a random num- 
ber normally distributed. The interpolation of V is done 
by means of a cubic spline connecting M random num- 
bers. Two extra conditions are needed in order to define 
a spline: we have taken F{i, 0) = and M) = 0. In 
our simulations M = 100. When the line advances be- 
yond Ui — M, a new spline, with M new random numbers 
is generated. Potentials for different i are independent. 

The value w is the center of mass of a confining po- 
tential for each point i, of the form '^{w — Ui)'^. In the 
simulation, w is increased from 0. For each value of w a 
metastable state is computed. Increasing w, a stationary 
sequence of metastable states (independent of the initial 
configuration) is reached, as observed in Ref. |48]. This 
is the steady state on which we focus. Our main results 
concern an elastic string in d = 1 of size L with RF dis- 
order, but we have also studied RB disorder, see Fig. [T| 



and a 2-dimensional elastic interface of size with peri- 
odic boundary conditions. As expected, for the depinning 
transition, the RB case falls in the same universality class 
as the RF case and results are very similar. 



III. DEFINITIONS AND OBSERVABLES 

For given w — wq the manifold moves to a metastable 
state Utop(a;), i.e. a state dynamically stable to infinites- 
imally small deformations. Following the notation of 
[ini [SD] , we define the center of mass of the metastable 
configuration 



u{wo) ■■= 



1 

I? 



dx 'i 



(3) 



with L the linear size of the system (number of points), 
and d the dimension. One then increases w, and a smooth 
forward deformation of u^{x) results (for smooth short- 
scale disorder) while the state remains stable. At some 
w — wi the state becomes unstable and the manifold, 
for w = Wi moves until it is blocked again in a new 
metastable state Uw^ix) (also locally stable). This pro- 
cess is called an avalanche and its size S is defined as 
the area swept by the line as it jumps between the two 
consecutive metastable states: 



S -.^ L"^ [u{wi) ~ u{wo)] 



(4) 



The distribution of avalanche sizes is expected to exhibit 
universality, i.e. independence of short scales, for sizes 
•5 > Snun- The short-scale cutoff Smm corresponds to 
the area spanned by a single monomer on the scale of 
the discretization of the disorder (in our units Smin — 1). 
In the limit m = a critical point is reached, resulting in 
a power-law distribution of avalanche sizes. To properly 
define the problem, including the stationary measure, it 
is essential to consider a small m > 0. Then, the internal 
correlation length is finite: it can e.g. be measured 
from the structure factor leading to |15] 



5/to 



(5) 



i„i is large in the small-m regime considered here. As 
a result, the distribution of avalanche sizes is cut off by 
the large scale Sm ^ S'min, defined as 



2(^) 



(6) 



It is expected to scale as Sm ^ Lf^'^ ^ m^'^^^ at small 
m. Here and below we define the (normalized) distribu- 
tion of avalanche sizes P{S), as well as its moments 



(7) 



from the sequence of measured avalanches Si, i 
l,...,iV. 
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FIG. 1: Random Field and Random Bond (RB) (d = 1). A 
fit witli a power law gives the exponent t — 1.08 ± 0.02. The 
agreement with Eq. ( 12 I is discussed in the text. 



The scale S„i is important as it allows to define univer- 
sal functions. In the variable s := S/S„i the avalanche- 
size distribution should become universal. Indeed, one of 
the predictions of the FRG theory is that if the exponent 
T satisfies 2 > r > 1 which is the case here, then the dis- 
tribution of avalanche sizes for S 3> 5'min takes the form 
as m ^ 0, i.e. 5,„ > S'min, 



P{S)dS := —p [-E-j 



(8) 



The function p{s) is universal and depends only on the 
space dimension d. Note that the normalized probability 
P{S) depends on the cut-off S'min via the first moment 
{S) which cannot be predicted by the theory, hence is an 
input from the numerics. It is important to stress that 
while the function p{s) is universal and convenient for 
data analysis, it is not a probability distribution and is 
not normalized to unity. Rather, it satisfies from its defi- 
nition (|8| and using ([6| the two normalization conditions 



(s)p = j ds sp{s) = 1 

j dss^p{s)^2. 



(9) 
(10) 



Here and below we use the notation (s)^ to denote an 
integration over p(s) and distinguish it from a true ex- 
pectation value over P{S), denoted (...). 



IV. THE AVALANCHE-SIZE DISTRIBUTION 



The rescaled avalanche-size distribution can be written 



as 



where r is the avalanche-size exponent, and /(s) the uni- 
versal cutoff function \57\ , which tends to a constant for 
s — > 0. For the present model, the only analytical pre- 
diction prior to our work |49l I50j concerns the exponent 
T, via the above mentioned NF [JTJ [55] conjecture 



' conj 



d+C ' 



(12) 



where C is the roughness exponent at the depinning tran- 
sition. Exact solution [5T] of a mean field toy model 
of avalanches, which turns out to be related to the fa- 
mous Gallon process |53| in genealogy, gives an exponent 
tmf = 3/2. This exponent is also the one expected if we 
replace d — due = 4 in the NF conjecture. This does 
however not constitute a first-principle calculation start- 
ing from the model of the pinned interface. The latter 
was only possible using the FRG [IHl \5U\ . The summa- 
tion of all tree diagrams within the FRG is shown to be 
asymptotically exact for d > 4 and leads to the mean- 
field prediction j49l [50] for r and for the full rescaled 
avalanche-size distribution (see below). 

We now discuss our numerical results starting with the 
avalanche-size exponent r. Note that the data in Fig. 
[T] contain both random-field and random-bond disorder 
and that, as expected from the universality of the depin- 
ning fixed point, the results are indistinguishable. Hence 
in the following we focus on RF disorder. For d = 1, a 
direct power-law fit of our numerical data (see Fig. [l]) 
gives 



= 1.08 ±0.02 



(13) 



This value has to be compared with the conjecture of 
Eq. (12 1. The roughness exponent is known numerically 
with a good accuracy from system sizes {L ~ lO'^) and 
TO = 0, as C = 1-26 ± 0.01 SD]. This value for C gives 



1.115 ±0.005. Hence the estimate (13 1 is slightly 



smaller than the value of r obtained from the conjecture. 
There are several possible explanations for this. 
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FIG. 2: Numerical extrapolation of the exponent to mass 
m = 0. We find C = 1-19 ± 0.01 
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First one notes that although (13) is extracted from 
pieces of p(s) which have already well converged in terms 
of TO and L, the resulting window of sizes is limited. Al- 
though we took this into account in estimating (13), we 



cannot exclude a further small upward shift in the central 
value as the window size increases. 

Second, we have also measured the effective C, exponent 
for the sizes and masses used here. From measurements 
of Sm we extract Q = 1.19±0.01 as can be seen on Fig.|2] 
We have checked that comparable estimates can be ex- 
tracted from the structure factor S{q), as measured also 
in [IS], using fits taking into account the mass. Inserting 
this value for an effective C, into Eq. ( [l2| ), this results in 
an effective value for Tconj — 1.086 ± 0.004, which is in 



much better agreement with our measured value (13). 

Finally, deviations from the conjecture for the asymp- 
totic value of T are still, strictly speaking, possible, but 
if they exist they must be around or below the error of 
0.02 in Eq. (13). This does not rule them out since, as 
discussed in [52], if present they are expected to be small 
[58|. 

Within the FRG UHl ISU] it is possible to compute the 
universal scaling function. For d > 4 summation of all 
tree diagrams gives 



/mf(s) 



1 



s/4 



The one-loop FRG calculation gives 



with exponents 



2^ 



exp ( C ^fs - 



1 



B 



(i-Ci)e 



(14) 
(15) 

(16) 
(17) 



where a = — 1(1 — Ci)^ and = 1/3 for the RF class, 
relevant to the present study. The constants A, B and C 
depend on e, and must satisfy the normalization condi- 
tions ([9]), (10). At first order in e they are C = — i^/7ra, 
B = l-a{l + ^), A = l + \{2~ 37E)a, 7e = 0.577216. 
As usual, the one-loop results for the exponents t, 5 and 
for the parameters A, B, and C are exact up to O(e^). 

To analyze our numerical data for the avalanche-size 
distribution, we have first computed (5*^) and {S) from 
the data, which allowed to determine numerically the uni- 
versal (and parameter-free) function p{s) using (|6| and 
([s]). Hence by construction the numerical data satisfy 
conditions ^ and (10). They are plotted in Figs. |3] and 
|4] with emphasis either on the power-law region or on 
the tail. Note that for the different values of m and L 
used here, the data have converged, with the exception of 
the last point for very large avalanches over-suppressed 
by the finite size of the interface in the smallest samples, 
and of the region of very small avalanches, which are cut 
off at s « l/Sm- 
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FIG. 3: Random Field {d = 1). Blow up of the power- 
law region. The red solid curve is given by Eq. ( 14 1 , the 



black dashed Une by Eq.([l5|, with A = 0.852, B = 1.56 and 
C = 0.56. 
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FIG. 4: Random Field (d = 1). Blow up of the tail region. 
The red solid curve is given by Eq.( 14 1 , the black dashed line 



by Eq.( 15 1, with A = 0.852, B = 1.56 and C = 0.56. 



To compare the numerical data with the mean-field 
and one-loop predictions, we use two procedures: 

In the first procedure we compare directly the cut-off 
functions /(s), see Fig s. [3] and |4) They are defined as 
/(s) :— s^p(s) where r is respectively T^uni — 1-08 for the 
numerical data, tmf = 3/2 for the mean- field prediction, 
and Tpade = 5/4 for the simplest Fade approximant of 
the one-loop result, i.e. setting e = 3 in ( [T6| . For d — 1, 
due to the large value of e = 3, the function Pi-ioopis) 
with T = 5/4, S = 7/6, A = 5/6- 7e/4, 5 = 5/3-1- 7e/6 
and C = y^/S does not have the correct normalization. 
We chose to introduce two rescaling factors 



p{s) = CiPi_ioop(c2s) 



(18) 



in order to enforce the conditions ([9| and (10). This 
procedure only changes the values of A, B and C in a 
consistent manner, see Figs. |3] and |4] Note that even 
though only /(s) is plotted, the chosen value of r changes 
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FIG. 5: Random Field {d = 1). Blow up of the power- 
law region. The red solid curve is given by Eq.(19l , the 



black dashed hne by Eq.^lSf, with A = 0.947, B = 1.871 and 
C = 0.606. 



FIG. 7: Random Field (d = 2, L = 100 for m = 0.1, 0.05; 
L — 200 for m = 0.025). A fit with a power law gives access 
to the exponent r = 1.3±0.01. The comparison with Eq.l 12 1 
is discussed in the text. 
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FIG. 6: Random Field {d = 1). Blow up of the tail region. 
The red solid curve is given by Eq.( 19 1 , the black dashed line 
by Eq.(fT5l), with A = 0.947, B = 1.871 and C = 0.606. 



the values of A, B, C via the normahzation conditions, 
hence must be discussed accordingly. 

A second approach, shown in Figs. |5] and [6] consists in 
fitting the same numerical curves as in Figs. |3] and |4] with 
either (i) an exponential function ("exponential fit") or 
(ii) the one- loop function ("fit one loop"), but using the 
numerically obtained exponent t = r^um ~ 1.08. The 
exponential fit reads 



(1 - - 



r(2-r) 



exp -1 



(19) 



All coefficients are determined as a function of r by the 
normalization conditions ([9]) and (10 1. Note that this 
exponential fit is mostly a guide to emphasize the sub- 
exponential tail apparent in the data. Similarly, for the 
one-loop fit we adopt the procedure described in the pre- 
vious paragraph, with r = Tnum everywhere instead of 



the one- loop Pade value r = 5/4. We expect this fit to 
be less sensitive to the lack of precision in the one-loop 
estimate of r for the large value of e = 3 relevant here, 
and to better capture the tail region. This is indeed what 
is found, see Fig. [6] It confirms the sub-exponential tail 
exponent (5 ~ 7/6 to a rather good precision. We stress 
that our procedure is not a fit using A,i3,C as fit pa- 
rameters, but that all parameters are specified by the 
one-loop prediction. 

We now turn to a 2-dimensional interface. The univer- 
sal function p(s) for c? = 2 and RF disorder is plotted on 
Fig.|7] From a direct power-law fit, we find 



_d=2 



1.3 ±0.02 



(20) 



This value has to be compared with the conjecture of 
Eq. ( 12 1. The roughness exponent at the depinning tran- 



sition is known numerically as Cnum — 0.753 ± 0.002 
which gives Tf=^^ = 1.2735 ± 0.0005. Although our value 
(20 1 of T is compatible with the conjecture, the precision 



is insufficient to conclude on possible small deviations 
from the latter. The mean-field and one-loop predic- 
tions discussed above are plotted for comparison, using 
the simplest one-loop Fade approximant, i.e. a = —4/9, 
T = 4/3, (5 = 10/9. After the above described procedure 
( 18 1 using the normalization conditions (|9]) and ( 10 ) this 
led to the values A = 0.92, B = 1.416 and C = 0.383. 



V. UNIVERSAL MOMENT RATIOS r„ 

Important universal quantities characterizing the 
avalanche statistics are the following universal ratios of 
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FIG. 8: Numerical values of the ratios r„ for bare random- 
field disorder {d = 1, L = 2000), and m = 0.02, m = 0.01. For 
m = a pole is expected for n — t. For n < t the divergence 
manifests itself through the data's dependence on m . 
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FIG. 10; Random Field {d = 1, L = 2000). Moment ratios 
r,i: blow up of the tail behavior. The mean field behavior is 
given by Eq.(|22|. The Fade (1,0) by Eq([23| and the Fade 
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FIG. 11: Random Field (d = 1, L = 2000). Moment ratios 
r„; blow up around the pole and comparison with Eq. (24 1. 



avalanche-size moments: 



(21) 



Here n can be non-integer. As shown in |50j all non- 
universal scales disappear in the ratios r„. Our numeri- 
cal findings are summarized in Fig. [8] The pole expected 
at n = T w 1.08 in the limit of infinite Sm/Smin man- 
ifests itself in a non-convergence of the numerical data 
upon lowering m. This is an independent method for 
calculating r. 

We now compare to the FRG calculation [5^. The 
function r„ can be evaluated in an e = 4 — cZ expansion. 
At the mean-field level (e — 0) 



2n- 1 
2n~3 ' 



(22) 



and a pole is found for n = tmf = 3/2. 

The one-loop e expansion leads to the following expres- 
sion l50l: 



r, 



l-loop 



3^ (2n-3)2r(n- I) 



(23) 

where Ci = 1/3 for RF. This expression corresponds 
to the Pade (1,0) in the e-expansion; we also use the 
Fade (0, 1). The comparison with the data is shown on 
Fig. [9] For the large-moment region, a blow-up is shown 
on Fig. [To] The agreement of the data with the two one- 
loop Pade approximants, as compared to mean field, is 
quite striking. 

However, both Fades break down close to n = r. We 
give another useful form for comparison to numerics. The 
idea is to isolate the simple pole which occurs in any 
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FIG. 12: The characteristic function Z{X} for d = 1 (left) and d = 2 (right). Mean field (soHd red/grey) from Eq. (291; one-loop 
from Eq. (301 (dashed black); and numerical results for L = 4000, m 
(291 and (IsOT for A : 



0.0005 (solid orange/green dots). The singularity in 
1/2 (indicated by a vertical dotted line) is smoothed out in the numerics, see figure 13 for details. 
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FIG. 13: Left: Log-log plot of dx In Z{X) versus A. For large A, we expect from Eq. ( 15 I that the analytic result (black, thick, 
dashed) obtained by integrating Eq. (281 numerically, has slope 1/(5 — 1) (blue dotted line is this asymptotics) . The numerical 
results are for L — 4000 and m = 0.02 (thick red) and m — 0.00125 (thin red). For large A, Znum(A) ~ ^ exp(A5'max), where 
Smax is the largest avalanche encountered in the simulation, and the curve saturates (dotted red lines). The larger Smax/Sm, 
the better the data. For m = 0.02 (thick red) this ratio is 16, whereas for m = 0.00125 (thin red) it is 10. The vertical dashed 
green line indicates the location A = 1/4 of the singularity in Z{X) at the mean-field level. Right: Slope of the function left, 
i.e. the effective exponent 1/(5 — 1). One sees that the effective exponent increases with increasing S'max/Sm. Our data, which 
clearly have not converged in terms of S'max/Sm, allow to estimate 1 < (5 < | from the maximum slope. 



dimension, as 



VI. THE CHARACTERISTIC FUNCTION Z{\) 



= \- Bn^d 



n — T 

Up to O(e^) corrections [50] 

l+TT . 



A, = 1 



Bn,d — 1 
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-(l-Ci)e 



(24) 



(25) 



7rr(n- i) - y/TrT{n-l) 



6(2n-3)r(n- i) 



^(l-Ci)e. (26) 



In Fig. [TTJthis formula is plotted setting e = 3. It shows 
that it works quite well, even close to the pole at n = r. 



It is useful to define a generating function of expo- 
nential moments, i.e. the characteristic function of the 
avalanche-size probability. Using the definitions (|6| and 
([t]), we define the normalized generating function Z{X) 



Z{\) 



s„ 



(S) 



XSi/S„ 



(27) 



By construction, Z{X) = A+A^+- • • . Since large negative 
A probe small avalanches, it is expected to be universal 
for A ^ — l/S'inin. In the universal range, its relation to 
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p{s) is 

poo 

Z{X) = / dsp(s) (e^" - 1) . (28) 
It has been calculated in |3H1 at the mean-field level: 



^mf(a) = ^ (i - VT^) 

At 1-loop order, it reads [121 [SU] 



(29) 



^i-ioop(A) = ^ (l - Vl - 4a) 



((3A + VI -4A - 1) log(l - 4A) - 2 (2A + Vl - 4A - l)) 
4V1 - 4A 



+ O {a') 



(30) 



where, as above. 



a = — "'"^^^ e = — |e. The comparison 
between theory and numerical data is presented on Fig. 
both for d 
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1, and d 



In these figures we have 
plotted (30 1, discarding the term 0{a^) and setting e = 3 



and e — 2 respectively. The plots show that the simplest 
extrapolation of the 1-loop correction is extremely good 
in calculating the behavior even for large negative A, as 
was already observed in the static case in [4^. It would 
be interesting to compare Z(X) for both cases numeri- 
cally. 

For large A, Z{X) is dominated by the largest avalanche 
S'max- If Smux/Sm ^ 1, then the tail-exponent S can in 
principle be extracted from the derivative of In Z(A), i.e. 
dx In Z{\) ~ A^/ ("5-1) in some window of A before it even- 
tually saturates to a constant ^ 5'max at larger A. Our 
data, which are plotted on Fig. [12] are not yet converged 
in terms of the ratio S'max/'S'm, but are sufficient to give 
the bound I < S < ^. One finally notes that although 
mean field works better for d — 2 than for d = 1, the 
1-loop corrections are necessary to account for the nu- 
merical data. 



VII. LOCAL AVALANCHE-SIZE 
DISTRIBUTION 



m in the definition of (f) (see above) such that S',„ and S!j^ 
both scale as m~'^~'^. 

For a more convenient comparison with numerics, we 
adopt a slightly different normalization as in Ref. |50j . 
and chose to normalize using S*^ rather than Sm- We 
estimate below the ratio — S'^/ Sm which allows to 
go from one set of definitions to the other. Hence the 
(normalized) local avalanche-size distribution P'^{S'^) is 
expected to take the form 



{Sty 



(33) 



where the universal function P(j,{x) = a^^y^p'^{a^x) in 

terms of the one defined in [50] and called p"^ there. By 
construction p^{x) satisfies the normalizations (Jol and 



(10 1. Similarly we define the generating function as 



Z{\) 



which reads Z{\) 
fined in I5D1. 



(S''^) N 



1 ^ 



\Sf /St 



(34) 



a(jjZ'^{\l a^) in terms of the one de- 



In [50', we have considered the following definition of 
the size of a local avalanche Sa,: 



Here we also define: 



C0 

0.rr, . 



2(V 



(31) 



(32) 



Of particular interest is the cross-section with a co- 
dimension one hyper-plane i.e. 0(x) — 'm~^d{xi), or more 
generally, with a co-dimension d' subspace. This cross- 
section has dimension = d ~ d'. We have chosen the 
factor of m in the definition of such that Sm and S"^ 
both scale as m~'^~^ . For d = 1 we consider a point, i.e. 
d' = 1, drf, = 0. Note that we always chose the factor of 
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FIG. 14: Random Field {d = I, d^, = Q, L = 4000). From the 
fit we get = 0.39 ± 0.01. 
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At present time we have only three analytical re- 
sults available to compare the numerical data on local 
avalanche-size distributions. First the conjecture put for- 
ward in [5D] and which generalizes ( 12 1 reads: 



(35) 



where we recall that = d — d! . Our numerical data for 
the point on a d = 1 string (i.e. one monomer) is shown 



on Fig. 14 and we find — 0.39±0.01. If we use the best 
present estimate C = 1-26 ± 0.01 we find t™"-' = 0.413 ± 
0.01. If we use the value of C = 1.19±0.01 extracted from 
the scaling of Sm, we find T™"j = 0.32 ± 0.02. While 
the values are roughly consistent, the precision on C, is 
crucial for a precise comparison. Inverting the conjecture 
(351, the measurement of — 0.39 ± 0.01 leads to a 
conjectured ^ = 1.24 ± 0.01. For a detailed discussion of 
the possible artifacts we refer to the discussion in Section 

El 

The second result is the exact expression of Z{\) and 
P(j,{s) in mean field, i.e. for d > A and d' = I. This in- 
volves a non-trivial summation of momentum-dependent 
tree diagrams using instanton calculus. It yields [50 that 
Z{X) is given by the solution of 



(Z-3)Z(2Z-3) ^ 9A 



(36) 



which vanishes at A = 0. This yields the series expansion 
Z{X) = A + A2 + f A3 + f A^ + ^A5 + 0(A6). We have 
compared this mean-field prediction and the numerical 
results in d = 1, (j^ = on figure 15 It is clear that loop 
corrections, yet to be computed, will play an important 
role, as was the case for bulk avalanches, see Fig. 
The function p^{s), as defined here, is found to be 
in mean field (i.e. at tree level) 



Ki 



(2V3) 



27rs 



(37) 
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FIG. 16: l/a^ as defined in equation ( |38[ ) and discussed below. 
For m <C 5/1/, ~ Lm. 



It satisfies the normalizations ^ and ([To| , and is related 
to Z(A) via Z{\) = /g°°dsp0(s) (e^'* - l). 

Finally, the mean-field calculation [5(1 also gives: 



Sm. 



1 

4 



(38) 



Corrections at 1-loop order slightly decrease this ratio. 
The e expansion predicts [54j 



1 

4 



a 1 



7tt 
6VS 



(39) 



with a — — f (1 — Ci) a — —2/3 here. Using the two 
Fade approximants gives the estimate l/a^ = 3.74 ±0.01 
which is consistent with the numerically observed value 
of l/aT'" = 3.7 ± 0.05 for L = 4000, m = 0.00125. 



VIII. CONCLUSION 



z 



0.5 


I 


-10 


-8 


-6 


-4 


-2 


0^0.5 














-1.0 














-1.5 














-2.0 





X 



FIG. 15: Z{\) both at the tree-level (solid/red), and numer- 
ically (green/orange dots), for RF disorder, m = 0.00125, 
d = 1, d,^ = 0. 



In this article, we have compared the numerically ob- 
tained avalanche-size statistics at the depinning transi- 
tion with the recent predictions from the functional RG 
based on an e = 4 — d expansion. The critical point of 
the depinning transition for an interface of internal di- 
mensions d — I and d — 2, driven quasi-statically in a 
random landscape in presence of an external quadratic 
well of curvature m^, is reached in the limit of to ^ 0. 
We have shown that the avalanche-size distribution P{S) 
takes the expected scaling form with the upper cutoff 
scale Sm ~ TO~''~^^, involving a universal function p{s) 
in the rescaled variable s = S/ Sm- As we confirmed, it 
does not depend on whether the microscopic disorder is of 
random-field or random-bond type. We have computed 
numerically the function p(s), its moments and its char- 
acteristic function and found in all cases good to excellent 
agreement with the predictions of the 1-loop FRG based 
on the extrapolation to d = 1 and d = 2. We have also 
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studied, for d = 1, the ZocaZ avalanches, and there too, we 
found a rather satisfactory agreement with available an- 
alytical predictions. However, it remains an outstanding 
challenge to compute the local avalanche-size distribution 
within the FRG beyond mean-field. 

Some fine points deserve discussion and further study. 
First we have not found any clear-cut signature that the 
conjecture for the avalanche exponent r be violated at de- 
pinning. However, we can not rule out such a violation 
below a 0.02 precision in r. A better numerical deter- 
mination of T, comparable in precision to the one which 
exists for the roughness ^ at m = would be crucial to 
confirm or invalidate the conjectured relation between r 
and Presently, the mass, i.e. the quadratic well, ap- 
pears necessary for a proper definition of the steady state, 
but unfortunately, this hampers the attempts at a more 
precise determination of t. 

Second, there has been a recent proposal, in the case 
of the random-field Ising model [55| , that avalanche-size 
distributions for statics and depinning are described by 
the same universal functions. Although the physics un- 
derlying this hypothesis is not clear to us, one may still 
ask the question for the present model [SS]. One may 
for instance compare our present results to the one in 
the statics in 49J. Currently, our precision is not suf- 
ficient to conclude. For instance, in c? = 2 and for RF 
disorder, the conjecture ([!]) for r gives r = 1.25 for the 
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tributions in each class. 



